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Figure 1 



Figure 2 



Figure 3 



Figure 4 



Station plot of MIZPAC 77 cruise. The position 
of the ice margin at the time of observation is 
also shown. 

Distribution and intensity of f i nestructure during 5 

MIZPAC 77. Symbols are described in the text. 

Isotherms (°C) are the maximum temperature in the 
water column. 

Station 80 salinity and temperature. The solid 10 

line was obtained while lowering, the dashed line 
while hoisting. The dotted line is the result of 
despiking the down trace but not smoothing. See 
Figure *4 for the smoothed curve. 

Station 80. The down-going salinity profile is 14 

on the right as a solid curve. The dashed curve 
is the result of despiking and smoothing. The two 
curves to the left give the contributions of the 
short time constant and long time constant correc- 
tions to the despiking. See Figure 3 for the un- 
smoothed corrected curve and to see the temperature 
transients which caused the salinity spikes. 



THE OCEANOGRAPHIC CRUISE OF USCGC BURTON ISLAND 
TO THE MARGINAL SEA- ICE ZONE OF THE CHUKCHI SEA- 

MI ZPAC 77 

by 

Robert G. Paquette and Robert H. Bourke 
Naval Postgraduate School, Monterey, CA 



I. INTRODUCTION 

This report presents the data and briefly describes the oceanographic 
results of the cruise of USCGC BURTON ISLAND into the region of the sea- 
ice margin of the Chukchi Sea during the period 2b July to 6 August 1977 
as part of the program designated MIZPAC 77. The primary objective of the 
cruise was to find and characterize f i nestructure in the vertical tempera- 
ture profiles and to discover its horizontal distribution and causes. This 
is the fifth cruise devoted to this general problem. Other cruises in 1971, 
1972, and 197^ were reported by Paquette and Bourke (1973, 1976) and MIZPAC 
75 by Zuberbuhler and Roeder (1976). 

I I . GENERAL DESCRIPTION 

The scientific group boarded BURTON ISLAND at Nome, Alaska on the 
afternoon of 2b July, one day later than in initial planning. The scientists 
and their affiliations were: 

Dr. Robert G. Paquette, Naval Postgraduate School, Chief Scientist 

Dr. Robert H. Bourke, Naval Postgraduate School 

LCDR Gordon P. Graham, Canadian Forces, Student at Naval Postgraduate 
School 

Mr. Jonathan D. Trent, Naval Postgraduate School 

Mr. Peter Becker, Applied Physics Laboratory, University of Washington, 
Seattle. 

The scientific group disembarked at NARL, Barrow, Alaska on the afternoon 
of 6 August. 

The measurements made were salinity and temperature profiles through- 
out the entire water column at 157 stations, using the Applied Physics Labora- 
tory portable, hand lowered CTD. The lowering rate of the CTD was about 1 m 
sec resulting in a data rate of approximately three points per meter. The 
latter was checked systematically with Nansen bottles lowered on a second 
wire. Prior to leaving each station, the temperature and salinity were 
plotted utilizing a Hewlett-Packard 9100 series computer/plotter system. 

These rough plots were used to make immediate assessments of the presence 
of fi nestructure and to aid in the decision of where to make the next few 
stations. 



A current meter, which was intended to orient lines of closely spaced 
stations along the flow direction, was found to be of little utility be- 
cause of the effects of the ship's iron on the magnetic compass. However, 
the failure to orient sampling lines along a presumed direction of propa- 
gation was at least partially overcome by running east-west lines as well 
as north-south lines in areas containing f i nes tructure . 

Because of the loss of one day of cruise time, no measurements were 
made in Bering Strait and the southern Chukchi Sea until the latitude of 
Pt. Hope was reached. Otherwise, we proceeded as in the general plan, ex- 
ploring for fi nes tructure along the ice margin from 70°N toward Barrow and 
studying it both in longitudinal and lateral sections when found. This re- 
sulted in much intensive study near Barrow followed, near the end of the 
cruise, by exploration near 71 °N and, finally, a few more measurements near 
Barrow. 

III. DATA 

The CTD was standardized by means of a Nansen bottle lowered on a 
second wire to a depth just above the sea floor. Forty-nine such compara- 
sons were in sufficiently unchanging water for temperature standardization 
and kO for salinity. The mean CTD temperature and salinity errors were 
0.08°C and -0.021 °/ 00 , respectively. The standard errors of these means 
were 0.0051°C and 0.005°/oo and the standard deviations were 0.036°C and 
0.025°/oo. A correction was applied to the depth data to account for the 
difference in density between fresh and salt water, as the CTD pressure 
sensor was calibrated in fresh water. 

The CTD records its data on a cassette which is eventually transferred 
to a seven-track tape by APL-UW for data editing and analysis at NPS. A 
computerized editing routine was written to remove erroneous data, inter- 
polate data where necessary, make temperature, salinity and depth corrections, 
and remove spurious salinity spikes. The despiking and editing procedures 
are described in some detail in Appendix A. 

Heading data for each station are listed in Appendix C. These contain 
station position and number, date/time of CTD lowering, water depth, type 
of navigation, wind, wave, and air temperature data, etc. Appendix B is 
an explanation of the codes used in Appendix C. 

Plotting routines were used to display property profiles for each 
station: temperature, salinity, sound speed, and density (cr^). These are 
compactly plotted four stations per page and displayed in Appendix D. Sta- 
tions taken in the deep water of the Barrow Canyon are shown two per page. 

IV. RESULTS 

The array of stations occupied is shown in Figure 1 together with an 
ice-margin position based principally upon observations made at the times 
stations were occupied. The ice-margin is thus not a single synoptic view, 
but a progressively distorted one which is more useful in describing ice- 
related phenomena. Synoptic views are also available. 
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Figure 1. Station plot of MIZPAC 77 cruise. The position of the 
ice margin at the time of observation is also shown. 



An interesting feature of the ice in the region of Pt. Barrow is the 
large embayment to the WNW. Most of the satellite-derived observations 
show the embayment more or less closed near Pt. Franklin but the ship found 
an open passage. A projection of ice toward Pt. Franklin is a feature which 
persisted for several weeks. 

Two groups of measurements were made from a hovering helicopter, one 
extending about 40 nmi into the ice and another group of stations in more 
or less open water which supplemented the ship-based stations. Helicopter 
operations proved highly useful, but could not always be used where they 
were most desirable because of persistent poor visibility. 

The distribution of f inestructure is shown in Figure 2 together with 
isotherms of the maximum temperature in the water column. The intensity 
of the f i nest ructure, measured from the maximum peak-to-peak temperature 
fluctuation in the profile, is categorized as weak, moderate or strong. 

Weak structure, with amplitude 0.2° to 0.5°C, is indicated by a circle with 
a central dot; medium structure, with amplitude 0.5° to 1.0°C, is indicated 
by a circle with a cross; strong structure, with amplitude greater than 1°C, 
is indicated by a solid circle. Fi nestructure, so classified, is construed 
as being in or below the thermocline, not in the usually shallow upper layer. 
An open circle indicates fluctuations of microstructure size or the absence 
of any notable structure at all. 

We decided to include the "nose" classification this year, after the 
opposite decision in 1975, because of a feeling that it may be an inter- 
mediate step in forming some types of deeper fi nestructure. The phenomenon 
we have called a "nose" in the past, in its simplest form, is a near-surface 
temperature maximum presumably formed by the cooling of an upper, warm 
mixed layer by ice or the overriding of such a layer with cold, dilute ice- 
melt water from nearby ice. The noses this year differed from most past 
results in that they were sometimes as much as 30 meters thick instead of 
the customary 10-15 m. The nose may be complex, having various levels of 
fi nestructure within itself. This f i nestructure sometimes cuts into the 
usually sharp thermocline so as to make the depth of the thermocline doubt- 
ful and the nose can thus grade into a situation which we should classify 

as normal f inestructure. Further, it seemed illogical to insist that the 

cooling be only at the surface since an initially cooled surface skin could 
be subsequently heated by solar radiation or overriden by warm dilute water. 
Therefore, a nose is recognized 

a. if there is a well-defined surface layer 5 or more meters thick 
with a cooling at the surface of at least 0.5°C as compared to 
the maximum temperature in the layer. 

b. if there is a well-defined surface layer in which there is 
moderate or strong f inestructure. 

The nose is indicated in Figure 2 by a tab on either side of the circle. 

If there is no structure in the nose, the tab is open; if there is medium 

or strong structure in the nose, the tab is filled solid. 
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Figure 2. Distribution and intensity of f inestructure during MIZPAC 77. 

Symbols are described in the text. Isotherms (°C) are the 
maximum temperature in the water column. 



It will be noted that much of the f i nest ructure is concentrated be- 
tween the 2° and 4°C isotherms of maximum temperature in the water column, 
a conclusion reached in studying earlier data. However, this is far from 
being an absolute rule: there are regions with no structure between the 

2° and 4°C isotherms and one instance of strong structure near the 9° 
isotherm. This conclusion therefore, requires closer examination. 

The following general results are presently apparent. 

1. The ice margin crossings in the vicinity of 167°W were surprisingly 
like the corresponding measurements of MIZPAC 72. This appears to be a 
region with little fi nest ructure. 

2. Much f i nestructure was found in an embayment in the ice roughly 
AO nmi square WNW of Barrow, even under ice. This is in agreement with 
the findings of MIZPAC 71 which operated in much the same area, but with a 
poorer understanding of the ice distribution. In the Barrow Canyon, where 
the water is deeper, fi nest ructure was found as deep as 50 m. This is 15 
to 20 m deeper than found on previous cruises. The fi nestructure and the 
embayment probably were due to the warm water of the coastal current which, 
in this case, had spread westward near Pt. Barrow. It is interesting that 
the ice was fairly close to shore farther south, near Pt. Franklin, and 

not melted away as one might expect if a continuous warm current moved north- 
ward near the surface. 

3. The coastal current was warmer than in the past, up to 8° to 10°C. 

It was close to the coast between Pt. Franklin and Pt. Barrow, but farther 
away farther south. This reinforces our feeling that the tendency of the 
current to be away from shore in the more southerly latitudes is a geo- 
graphical rather than a temporal phenomenon. Five crossings of the coastal 
current, between Icy Cape and Pt. Barrow, were made during this cruise. 
Analysis of these crossings should provide insight into the route of the 
warm Alaskan Coastal Current as it flows northeastward from the shallow 
waters off Icy Cape to the deeper waters of Barrow Canyon. 

4. The salinity of bottom water behind the ice edge was relatively 
high, 33.^"33.6°/ 00 , like in MIZPAC 71 and MIZPAC 72. What causes this 
large year-to-year variation in mean bottom water salinity, often greater 
than 0.5°/oo»has not been established. It could be controlled by the amount 
of brine rejection occurring in winter in the shallows of the Chukchi Sea, 

a phenomenon which is a function of the severity of the winter. However, 
the cause might also be in river flows into the Bering Sea during the pre- 
vious summer. 

At the extreme bottom of many casts it was noted there was an apparent 
increase in salinity in the last half meter above bottom. We are undecided 
about whether this increase is an instrument artifact or a real phenomenon. 
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APPENDIX A 

DESP IKING AND DATA EDITING 



Introduction 



Salinity spikes appear to be considerably more serious in 1977 than 
they were in 1975. This must be due to the instrument modifications made 
in the interim. Some objective despiking routine obviously was required. 

In October 197** we had tested a despiking program for an old model Plessey 
STD based entirely on the concept of a first-order response equation for 
temperature and some phase shift between salinity and temperature, which 
Scarlet reported using, without the phase sh i ft , in 1975. The method did 
not work well on our data. Constants which would correct the top half of 
a spike would overcorrect the bottom half. The reasons may lie in second- 
order effects, perhaps complicated by the fact that we worked from the 
graphically recorded output rather than from digital data and thus included 
recorder response in the data. Or it may be that the technique which Scarlet 
found suitable for small spikes could not handle the very large spikes 
occurring in the Chukchi Sea data which were further exaggerated by the long 

thermometer time constant of the STD. Nevertheless, necessity prompted us 

to try again, with some success. 

The data editing program to be described also is not completely new. 

We had one in MIZPAC 75. But the present program is considerably more 

versatile and automatic. It contains a "ratchet" or "latch" subroutine 

which prevents depth reversals similar to the one described by Scarlet. 

Our routine differs in that Scarlet threw out offending data whereas we 
replace offending points with interpolated points. This is necessary be- 
cause we do not record time on the tape and must assume that the time step 
is uniform when despiking. 

Despiking experiments have proceeded at the Applied Physics Laboratory, 
University of Washington (APL) and at the Naval Postgraduate School (NPS) 
with considerable general communication between the two laboratories. The 
method of handling the short time constant below is due mainly to G. 

Garrison of APL, the method of correcting temperature and the long time 
constant is due to R. G. Paquette of NPS. The details of implementation 
by APL and NPS are different and we have not i ntercompared results. The 
editing program described was evolved at NPS. APL uses a different technique. 

Facts about the Instrument 

The CTD is a portable digital instrument with a conventional three- 
electrode conductivity cell and a thermistor as a temperature sensor. The 
cell constant of the cell is remarkably stable, judged by classical ideas 
about polarization phenomena in cells with unplatinized electrodes and 
relatively high electrode current densities. Perhaps this is due to the 
high operation frequency, ca_ 10 kHz, and the low electrode voltage, 13 mV. 

The three sensors of the instrument produce electrical frequencies which 
are counted in the instrument and the counts are recorded digitally on 
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cassette tape. The sensor data set is sampled about 3 times per second 
in the order conductivity, temperature, depth. The fact that the conduc- 
tivity and temperature measurements are separated by about 30 ms is a fac- 
tor in the generation of salinity spikes and their subsequent removal. 

Salinity Sp i kes- General 

Spurious salinity spikes in a CTD record result from the lack of 
simultaneity between the conductivity and temperature measurement. The 
spikes are of consequence only when there are sharp temperature transients. 

A time lag may occur because of sampling lag, as above, or it may occur be- 
cause of response lags in the sensors. The effect exists whether the salinity 
is computed digitally or by an analog circuit. When the temperature is 
retarded compared to the conductivity by either mechanism, the spike has 
the opposite sense from the first derivative of the temperature- t i me curve 
and conversely. If both measurements are effectively simultaneous but both 
are wrong because of sensor lag, the salinity, if slowly changing, may be 
correct . 

Sources of spikes 

In our investigation of the cause of salinity spikes, which sometimes 
exceeded l°/ 0 o in magnitude, we found several sensor response errors. 

° A first-order lag in the temperature, with time constant 0.05 sec. 

° A first-order lag in the flushing of the conductivity cell with length 
constant about 18 cm. 

° An error due to non-s i mu 1 taneous sampling of conductivity and temperature 
and to a physical vertical displacement between the conductivity cell 
and the temperature sensor. 

° A long time constant in conductivity response apparently due to heat 
storage somewhere in the cell structure. 

Recognition and Evaluation of the Long Time Constant 

It may be of interest to know how the long time constant came to be 
recognized and evaluated. It was found only because data were recorded on 
the upward traverse of the CTD as well as when it was going down. It was 
then noted that the down-going salinity curve often did not agree with the 
up-going curve in the region below the thermocline where the temperature 
was unusually constant. No temperature- i nduced salinity error should have 
been generated in situ . Obviously the error must have been a long-lasting 
effect of a salinity or temperature transient which occurred at a shallower 
level and affected the down trace. The upward-going trace must have been 
essentially correct (except for the neglected effect of salinity change on 
conductivity) so long as the temperature was constant because the sensors 
usually were held near bottom for perhaps 10 seconds before beginning to 
hoist. The phenomenon under discussion may be seen in Figure 3, which shows 
the down and up salinity traces in Station 80 as well as the down and up 
temperature traces and the corrected salinity trace. One may easily find a 
time constant of about 5 sec in the approach of the salinity down trace to 
the up trace below the thermocline, taking the lowering rate to be 1 m/sec. 
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Figure 3. Station 80 salinity and temperature. The solid line was obtained 
while lowering, the dashed line while hoisting. The dotted line is the 
result of despiking the down trace but not smoothing. See Figure 4 for the 
smoothed curve. 



If such a time constant exists in other STD's and CTD 1 s , there must be 
difficulties in removing spikes by the use of a single time constant. This 
is because the short and long time constants contribute differently to narrow 
spikes than to broad spikes; the short time constant influences the former 
more than it does the latter and conversely for the long time constant. 

The Possibility of Second-Order Effects 



All of the discussion of spike correction below will assume the simple 
first-order response equation for sensors. It should be kept in mind, 
however, that substantial second-order components may be present. Whenever 
there are two thermal masses in series with two thermal resistances involved 
in the transfer of heat from the water to the interior of the sensor, a 
second-order term enters the response equation. The temperature and salinity 
sensors must have had these character i s t i cs to some extent, but it was not 
possible to find the multiplier of the second-order term. A second-order 
corrector basically takes out curvatures. It was found that the second-order 
constant could be adjusted to take out anomalies in the curves beautifully, 
but there was no way to tell whether the corrections corresponded to reality. 
Therefore, after a few experiments, the refinement was abandoned. It may 
be possible to revive it at some later date when more is known about the 
causes of variability of the sensor response constants. 

Lack of Constancy of the "Constants 11 



A major problem in devising a spike correction technique is the fact 
that the response constants of the sensors are not constant. They change 
with rate of lowering, wire angle and relative horizontal motion between the 
sensor head and the surrounding water. If the constants had been nicely 
fixed, there would have been s tra ightforward techniques for evaluating 
several constants. As it was, there was little hope of evaluating the 
second-order constant and little justification for trying to compute the 
effect of the flushing constant of the cell upon the conductivity directly, 
a natter which will be discussed below when the methods of spike correction 
are related in detail. 

Theory and Algorithms for Correction of Spikes - Thermometer Correction 

In the following paragraphs the nature, causes and method of correction 
of the several time constants are discussed. 



We assumed that the thermometer followed a simple first-order 1 aw 



T - T' 



k 



dT' 
T dt 



where T is the water temperature, T' the observed temperature, k^ the time 
constant of the thermistor and _t the time. The bare thermistor had a nominal 
time constant of 0.03 sec. The protective shroud was expected to add to this. 
After a number of trials we picked 0.05 sec as the largest value which 
would not produce an overshoot at the bottom of the sharpest thernocl i nes . 
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The temperature derivative was obtained from the first central difference 




where h_ is the time step between samples, 0.33 sec. 

The Time Lags and the Flushing Time Constant 

The difference in sampling times of temperature and conductivity was 
0.03 sec (temperature latest) plus an addendum of 0.01 sec because the 
thermometer was 1 cm deeper than the mouth of the cell. These were com- 
pensated by interpolating backward toward the previous temperature a frac- 
tional distance (0.01 + 0.03)/h or 0.012, using the equation 

AT. = k T. . + ( 1 -k ) T. - T'. 

J c j-1 c j j 

where the constant k would be expected to be 0.12, from above, but was 
more effectively about 0.65. The extra 0.53 can be shown to approximate 
a correction for the flushing lag in the conductivity cell, which apparently 
has a time constant of 0.53 h or 0 . 1 8 sec. This corresponds to a length 
constant of 18 cm, or one cell length, which is reasonable. 

There are some known theoretical weaknesses in treating the flushing 
time constant of the cell in this way. The method involves the assumption 
that there are no significantly large conductivity gradients that are 
not caused exclusively by temperature changes, which is the same as assuming 
that the true salinity always changes slowly. It also involves the assumption 
that the first backward temperature difference is proportional to 3c/9t, 
where c_ is the conductivity. These are not always good assumptions. However, 
reasonable success with the techniques led us to postpone further refinements 
to a later date. 

The Long Time Constant 



The long time constant in the cell was modeled successfully as a thermal 
mass coupled to true water temperature by a thermal resistance corresponding 
to a time constant k^. This mass passes fraction £ of the difference be- 
tween its temperature and true water temperature to the cell. It is a 
fairly realistic model. The thermal mass might be the cell body itself or 
the surrounding protective shroud. 



The temperature, 
s tart i ng wi th TC^ = T 
tween TC and T drives 



TC, of the thermal mass is accumulated at each step, 
and assuming that the mean temperature difference be- 
the temperature change in accordance with 



TC. = TC. . 

j j-i 



T ,+T . . - TC . +TC . . 

j j -1 j nL 



J k l 



which ultimately gives 
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